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I.  INTRODUCTION 


This  report  describes  a synthesis  procedure  for  efficiently  computing 
the  aperture  distribution  of  an  antenna  based  on  knowledge  of  its  measured 
radiation  pattern.  Although  the  methods  considered  are  generally  appli- 
cable to  the  problem  of  antenna  pattern  synthesis,  they  are  considered 
here  from  the  standpoint  of  providing  a useful  tool  in  the  computer-aided 
design  of  on-aircraft  antenna  systems.  As  a first  step  in  a typical 
on-aircraft  antenna  problem,  the  designer  is  confronted  with  the  task 
of  selecting  one  or  more  candidate  antenna  configurations  which  show 
promise  of  providing  the  desired  performance  in  the  on-aircraft  environment. 
The  selection  of  candidate  antennas  is  usually  made  by  analysis  of 
measured  pattern  data  taken  with  the  antenna  in  isolation  or  mounted  in 
a finite  ground  plane.  Traditionally,  the  antenna  designer  tests  each 
of  the  antennas  on  an  actual  or  scale-model  aircraft,  and  after  many 
measurements  a suitable  on-aircraft  antenna  system  eventually  evolves. 


In  recent  years  it  has  often  been  possible  to  use  computer  simulation 
to  carry  out  much  of  the  detailed  on-aircraft  antenna  design.  Such 
computer  aided  designs  have  typically  relied  on  assumed  approximate 
aperture  distributions  to  analytically  represent  the  candidate  antenna, 
since  the  detailed  distribution  for  the  actual  antenna  is  usually  not 
available.  The  synthesis  procedure  described  in  the  present  report 
provides  an  efficient  means  for  overcoming  this  difficulty. 


The  approaches  employed  are,  with  some  modifications,  an  adaptation 
of  basic  synthesis  techniques  presented  by  Mautz  and  Harrington  [1].  The 
viewpoint  here  differs  from  the  usual  synthesis  situation,  however,  in 
that  the  desired  pattern  is  assumed  to  be  actual  and  available  data,  rather 
than  an  idealized  specification.  More  specifically  the  procedure  as 
employed  here  requires  detailed  pattern  magnitude  data  for  the  candidate 
antenna,  taken  in  the  principal  planes  (i.e. , E-  and/or  H-planes)  of  the 
antenna.  Since  only  pattern  magnitude  data  is  assumed  to  Be  given,  an 
iterative  numerical  procedure  is  required  to  evolve  the  unspecified  pattern 
phase  information. 

The  methods  to  be  discussed  are  equally  applicable  to  continuous  or 
discretely  distributed  aperture  antennas,  and  the  antenna  aperture  may  be 
isolated  or  mounted  in  a finite  ground  plane.  For  simplicity,  the 
discussion  specifically  considers  the  antennas  to  be  linearly  polarized, 
with  a planar  aperture  mounted  in  a planar  ground  plane.  The  methods, 
however,  are  general  and  could  readily  be  extended  to  treat  arbitrarily 
polarized  antennas,  as  well  as  arbitrarily  curved  apertures  and  ground 
planes. 

The  following  section  presents  a brief  discussion  of  the  theoretical 
basis  for  the  synthesis  methods  employed.  Section  III  discusses  details 
of  its  numerical  implementation  and  presents  several  example  solutions 
obtained  for  specific  measured  antenna  geometries. 
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II.  THE  THEORETICAL  BASIS  OF  THE  SYNTHESIS  PROCEDURE 

This  section  discusses  the  general  theoretical  aspects  of  a 
synthesis  procedure  for  computing  the  aperture  amplitude  and  phase  distri- 
bution of  an  antenna,  given  its  far-field  magnitude  pattern.  The 
theoretical  approaches  employed  are  similar  to  those  reported  by  Mautz 
and  Harrington  [1].  The  subsequently  discussed  (Section  III)  numerical 
implementation  of  the  basic  procedure  is  more  specifically  geared  to  our 
application  requirements  as  outlined  in  the  Introduction. 

The  general  problem  to  be  solved  is  depicted  in  Figure  la.  Starting 
with  a given  or  specified  radiation  pattern,  it  is  desired  to  find  an 
aperture  distribution  (in  both  amplitude  and  phase)  whose  field  suitably 
approximates  the  given  radiation  pattern.  The  aperture  distribution  may 
be  either  discrete  or  continuously  distributed  and  the  ground  plane 
containing  the  aperture  may  be  finite  (including  no  ground  plane)  or 
infinite.  For  the  purposes  of  analysis,  both  the  computed  source  (i.e., 
the  computed  aperture  distribution)  and  its  resulting  field  are  considered 
to  be  representable  by  discrete  sets  of  values.  This  situation  is 
depicted  in  Figure  lb.  If  the  true  source  is  continuously  distributed 
the  discrete  representation  is,  of  course,  an  approximation.  If  the 
true  source  is  discrete  the  discrete  representation  may  be  exact  in  theory 
(although  any  of  a multitude  of  factors  would,  in  practice,  prevent  an 
exact  realization).  The  N elemental  sources  (actually  their  relative 
weights)  representing  the  true  source  can  be  expressed  in  matrix  form  as 


f 


(1) 


In  a similar  manner  the  synthesized  pattern,  computed  at  M field  points, 
can  be  expressed  as 
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RADIATION  PATTERN 


SYNTHESIZED  — 
RADIATION  PATTERN 


A typical  antenna  geometry  and  its  analytical  description 

a)  Actual  antenna  geometry 

b)  Analytical  model  of  the  antenna 


Figure  1 


The  source  and  resulting  field  are  then  formally  related  according  to 
the  relation 


[T]f  = g. 


(3) 


Here  [T]  is  an  MxN  matrix,  expressed  as 


[T] 


= CTmn]  « 


|_[m1 tmn 

where 


T = the  field  of  the  nth  elemental  source,  evaluated 
at  the  mth  field  point  location. 


(4) 


The  synthesis  problem  can  now  be  expressed  by  the  approximate 
relation 


[T]f  % g0  (5) 

where  g^,  is  the  finite  or  infinite  set  of  specified  field  values  (Note 
that  g in  Equation  (3)  is  the  set  of  M synthesized  field  values,  resulting 
from  the  N elemental  sources).  It  may  be  noted  that  there  are  several 
ways  to  solve  Equation  (5)  for  the  unknown  set  of  sources  [fn].  For 
example,  if  the  matrix  [T]  is  square  and  non  singular  Equation  (5)  can 
be  inverted  directed  to  yield  the  set  of  values  [fn].  In  this  case  the 
field  can  be  synthesized  correctly  at  M pattern  points,  but  no  control 
of  the  pattern  is  obtained  between  these  points.  If  M is  less  than  N, 
the  problem  is  underspecified  (i.e.,  there  are  more  unknowns  than 
independent  equations),  and  it  is  possible  to  obtain  more  than  one  set 
of  [fn]  which  satisfy  Equation  (5).  Finally,  if  M is  greater  than  N, 
then  an  exact  solution  usually  doesn't  exist,  but  a unique  least  squares 
solution  can  be  found.  This  latter  solution  case  is  usually  the  most 
common  one  encountered  in  practice  and  is,  therefore,  the  one  adopted 
in  this  study. 

It  can  readily  be  shown  (e.g.,  see  Reference  [1])  that  the  set  of 
[fn]  can  be  determined  which  minimizes  the  pattern  error  in  the  sense 
of  obtaining  the  least  mean  square  error  when  Equation  (5)  is  solved 
according  to  the  relation 
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(6) 


[f„]  - {tfHTlj-'tf]  g0. 


where  the  tilda  over  the  matrix  indicates  conjugate  transpose.  The 
mean  squared  pattern  error  resulting  from  the  above  solution  is  given  by 


M 

e=  I 

m=l 


N 

y f-  T - h.  e 


J6 


m 


n=l  " m 
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(7) 


where  the  mth  specified  pattern  value  has  been  expressed  as 


s • 


’om 


(8) 


The  error  as  expressed  by  Equation  (7)  assumes  that  the  pattern  is 
specified  in  both  amplitude  and  phase  (Equation  (8)).  As  discussed  in 
the  Introduction,  pattern  phase  data  may  not  be  available  or,  for  other 
reasons,  mav  not  be  specified.  In  this  situation  a solution  to 
Equation  (6)  which  minimizes  the  pattern  error  can  still  be  obtained. 
Specifically,  for  any  set  of  solutions  [fn]  the  pattern  error  is 
minimized  when  the  terms  in  Equation  (7)  are  in  phase.  This  requires 
that  the  phases,  Bm,  of  each  of  the  field  values  satisfy  the  relation 


N 

L,  fn  ^mn 
n=  l 

"R 

l fn  L 

ril  n mn 


(9) 


When  pattern  magnitude  only  data  is  specified  (i.e.,  the  hm,  in 
Equation  (8))  an  iteration  process  is  used  to  find  a set  of  Sm  which 
satisfy  Equation  (9)  (and  which  thereby  yield  a solution  for  the  sources 
[fn]  which  minimize  the  pattern  error  e according  to  Equation  (7)). 

The  steps  in  the  basic  iteration  procedure,  as  given  in  Reference  [1], 
are  summarized  as  follows: 

1.  Assume  starting  values  for  the  field  point  phases  3], &2>‘  "&M- 

2.  Keep  the  8m  fixed  and  compute  the  elemental  source  values  fn 
which  minimize  the  pattern  error  e,  using  Equation  (6). 

3.  Keep  the  fn  fixed  and  compute  the  6m  which  minimize  e,  using 
Equation  (9). 

4.  Repeat  steps  2 and  3 in  iterative  fashion  until  a small  and 
stable  value  of  e is  obtained. 
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Although  numerical  computations  performed  with  this  basic  iteration 
procedure  lead  to  a convergent  solution  for  the  elemental  sources  fn, 
the  rate  of  convergence  has  been  found  to  be  rather  slow  in  some  cases. 

In  order  to  reduce  the  number  of  iterations  required  for  convergence,  a 
modified  procedure  was  developed.  Within  a given  iteration,  the  modified 
procedure  makes  incremental  changes  in  both  the  computed  fn  and  the 
computed  (i.e.,  following  steps  2 and  3,  respectively,  above),  based 
upon  the  change  in  these  values  between  successive  iterations.  The 
modified  procedure  is  included  in  the  computer  code  presented  in 
Appendix  A,  and  has  been  found  to  reduce  the  number  of  iterations  required 
by  about  a factor  of  ten  in  most  instances,  compared  to  the  basic  pro- 
cedure. Finally  it  should  be  noted  that  although  the  synthesis  procedure 
discussed  here  generally  yields  satisfactory  solutions,  the  solutions 
obtained  are  not  unique.  For  example,  depending  upon  the  assumed  starting 
values  for  the  field  point  phases  (i.e.,  step  1 in  the  iteration  procedure) 
one  could  obtain  somewhat  different  solutions  for  the  sources  fn. 


III.  NUMERICAL  IMPLEMENTATION  OF  THE  PROCEDURE 

This  section  contains  a discussion  of  some  of  the  specific  details 
pertinent  to  practical  application  of  the  synthesis  procedure.  Examples 
are  presented  which  demonstrate  the  techniques  and  are  representative  of 
results  obtainable.  For  the  present  discussion  it  is  assumed  that  the 
antenna  is  linearly  polarized,  and  that  its  pattern  is  specified  by 
measured  principal  plane  pattern  data.  Also,  the  size  and  shape  (e.g., 
rectangular)  of  the  aperture  are  assumed  to  be  known  and,  if  tne  antenna 
aperture  is  flush-mounted  in  a finite  ground  plane,  the  dimensions  of  the 
ground  plane  are  also  assumed  to  be  given. 

% 

With  the  antenna  characteristics  physically  and  electrically  defined 
as  indicated  above,  the  first  step  in  numerically  implementing  the 
synthesis  procedure  involves  selection  of  appropriate  values  for  M and  N, 
which  it  may  be  recalled  are  the  total  number  of  specified  field  point 
values  and  the  total  number  of  elemental  sources,  respectively.  In 
general,  appropriate  selection  of  the  values  for  both  M and  N are  directly 
related  to  the  size  of  the  antenna  aperture.  (That  is,  as  the  aperture 
size  increases,  the  appropriate  values  for  both  M and  N should  be 
correspondingly  increased.)  As  will  be  demonstrated  shortly,  adequate 
representation  of  the  aperture  distribution  usually  requires  that  the 
elemental  sources  be  spaced  on  the  order  of  .25A  apart.  With  the  size  of 
the  aperture  known  and  the  elemental  source  spacing  selected  as  just 
indicated,  a suitable  value  for  the  number  of  sources  N can  be  chosen. 

As  noted  earlier  the  number  M,  of  field  point  values  selected,  will 
generally  be  considerably  larger  (e.g.,  3-20  times  larger)  than  the  number 
of  elemental  sources  employed.  In  practice,  the  actual  number  of  field 
point  values  selected  involves  making  a tradeoff  between  desired 
accuracy  and  computational  efficiency.  On  the  one  hand  enough  field 
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point  values  should  be  selected  from  the  given  pattern  data  to  describe  the 
detailed  structure  of  the  pattern  with  reasonable  accuracy.  On  the  other 
hand,  selection  of  a large  number  of  field  points  increases  computation 
time  and  may  not  yield  a substantially  better  solution.  In  any  event, 
an  appropriate  selection  of  field  points  (both  in  number  and  relative 
location)  can  usually  be  obtained  by  inspection  from  the  given  pattern 
data. 

Having  selected  a list  of  field  point  values  and  an  appropriate 
number  of  sources  as  described  above,  an  analytical  model  of  the  radiating 
aperture  is  established  to  compute  the  elements  of  the  matrix  [Tmn] 
associated  with  the  fields  of  the  elemental  sources,  as  defined  in 
Equation  (4).  For  the  examples  to  be  presented,  the  elemental  sources 
employed  are  magnetic  line  sources.  A typical  geometry  depicting  their 
placement  within  the  known  aperture  configuration  is  depicted  in  Figure  2. 


Figure  2.  Geometry  depicting  computation  of  a matrix  element 
Tmn.  A matrix  element  is  the  computed  value  of 
the  total  field  radiated  by  the  nth  elemental 
source  fn,  evaluated  at  the  mth  specified  location. 
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As  may  be  recalled  from  the  previous  section  each  element  of  the  matrix 
[Tmn]  is  simply  the  computed  value  of  the  field  of  the  nth  elemental 
source  evaluated  at  the  mth  selected  field  point  location.  Evaluation  of 
the  Tmn  is  straight-forward,  and  need  be  computed  only  once  (assuming  that 
the  values  are  stored).  If,  as  depicted  in  Figure  2,  the  aperture  is 
mounted  in  a finite  ground  plane,  the  diffracted  fields  from  the  plate 
edges  may  be  included  in  the  computation  by  superposition  with  the  direct 
field  from  each  elemental  source.  When  the  effects  of  a finite  ground 
plane  are  included  in  this  manner  selection  of  field  point  values  will  be 
most  effective  if  all  of  the  field  values  are  chosen  from  the  half-space 
directly  illuminated  by  the  aperture  (i.e.,  the  half-space  on  the  antenna 
side  of  the  ground  plane,  in  Figure  2). 

As  a first  typical  example  of  the  type  of  antenna  which  can  be  treated 
with  the  foregoing  analysis,  the  antenna  configuration  depicted  in 
Figure  3a  was  constructed.  The  measured  E-plane  radiation  pattern  (i.e., 
pattern  plane  is  the  x-z  plane  in  Figure  3a)  for  this  antenna  is  depicted 
in  Figure  3b.  Since  the  1.33  inch  width  of  the  antenna  aperture 
corresponds  to  about  1.25  wavelengths  at  the  11  GHz  measurement  frequency, 
six  magnetic  line  sources  were  selected  to  model  the  aperture.  The 
sources  were  uniformly  spaced  and  oriented  parallel  to  the  y axis  (see 
coordinate  system  of  Figure  3a).  For  specified  pattern  data  points,  a 
set  of  91  field  values  were  selected  from  the  measured  pattern  data  of 
Figure  3b.  In  order  to  demonstrate  how  the  iterative  synthesis  pro- 
cedure converges  to  a final  solution  for  the  antenna  of  Figure  3, 

Figure  4 shows  a sequence  of  patterns  computed  after  1,  10,  20,  40  and  80 
iterations  of  the  numerical  synthesis  procedure.  The  measured  data 
points  employed  to  obtain  the  computed  solution  are  shown  on  each  pattern, 
and  the  computed  relative  magnitude  and  phase  of  the  six  elemental 
sources  (from  which  each  pattern  is  computed)  are  listed  adjacent  to  each 
of  the  respective  patterns.  Finally,  the  pattern  errors  obtained  at  each 
of  the  indicated  stages  of  the  solution  are  also  shown.  The  desired 
information  in  the  solution  as  depicted  sequentially  in  Figure  4 is,  of 
course,  the  list  of  complex  source  values  obtained  after  solution 
convergence  is  achieved.  Convergence  of  the  solution  is  achieved  when 
the  list  of  source  values  obtained  remains  relatively  constant  (or 
alternatively,  when  the  pattern  error  values  reach  a relatively  stable 
minimum  value).  For  the  solution  depicted  in  Figure  4,  a convergent 
solution  was  obtained  after  about  40  iterations  (Figure  4d).  This  can 
be  seen  by  comparing  the  results  obtained  after  40  and  80  iterations 
(Figures  4d  and  4e,  respectively). 

Although  one  might  choose  to  gauge  solution  convergence  visually  by 
comparing  the  synthesized  pattern  with  the  specified  measured  data 
points,  this  is  not  always  a reliable  approach.  Because  of  errors  in 
the  measured  pattern  data  (or,  equivalently,  failure  to  include  a minor 
feature  of  the  physical  problem  in  the  analytical  model),  the  synthesized 
pattern  may  exhibit  an  irreducible  error  in  some  region  of  the  pattern. 
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ANGLE  e ( DEGREES  ) 


Figure  3.  The  first  example  antenna  configuration 
and  its  measured  pattern. 

a)  The  antenna  geometry 

b)  Its  measured  E-plane  radiation  pattern 


MEASURED 
DATA  POINTS 


demonstrate,  in  sequence,  the  manner  in  which  a convergent  solution  is 
obtained  by  use  of  the  iterative  synthesis  procedure. 


(b)  Solution  obtained  after  10  iterations 


(c)  Solution  obtained  after  20  iterations 


obtained  after  40  iterations 


MEASURED 


(e)  Solution  obtained  after  80  iterations 


even  after  a convergent  solution  has  been  attained.  This  situation  is 
apparent  in  a portion  of  the  computed  patterns  of  Figures  4d  or  4e.  Since 
measurement  errors  tend  to  increase  as  the  pattern  levels  decrease,  it 
may  be  advantageous  to  specify  fewer  measured  data  points  in  the  lower 
radiation  level  regions  of  the  pattern.  This  technique  is  demonstrated 
in  Figure  5,  which  shows  the  solution  of  Figure  4 recomputed  using  fewer 
measured  data  points  (i.e.,  65  data  points  were  used  instead  of  the 
original  91).  In  selecting  the  measured  data  points  of  Figure  5,  the 
original  data  set  was  thinned  by  retaining  only  every  other  one  of  the 
original  points  in  regions  of  the  pattern  with  signal  levels  in  the 
range  between  10  and  15  dB  below  the  peak  pattern  level,  and  retaining 
only  every  third  data  point  in  regions  of  the  pattern  with  signal 
levels  more  than  15  dB  down.  It  may  be  noted,  by  comparing  the  solution 
of  Figure  5 with  that  of  Figure  4d,  that  virtually  the  same  solution  is 
obtained  when  fewer  data  points  are  employed  in  this  manner.  In  addition, 
the  computation  time  was  reduced  substantially  by  use  of  fewer  specified 
pattern  data  points. 

As  a second  application  of  the  synthesis  procedure,  the  antenna 
geometry  of  Figure  6a  was  considered.  This  antenna  was  constructed  and 
tested  several  years  ago  by  Walter  [2]  as  part  of  a study  of  endfire 
radiators.  We  have  applied  the  synthesis  procedure  to  an  available 
measured  pattern  for  this  antenna  as  shown  in  Figure  6b.  Although  the 
antenna  geometry  of  Figure  6a  is  similar  in  basic  design  to  that  just 
considered,  it  differs  from  the  latter  in  two  important  respects. 
Specifically  the  aperture  of  the  antenna  in  Figure  6a  is  roughly  three 
times  as  large  as  that  shown  in  Figure  3,  and  in  addition,  its  aperture 
and  waveguide  feed  were  filled  with  dielectric.  For  the  aperture  of  the 
antenna  of  Figure  6,  12  magnetic  line  sources  were  used  to  model  its 
aperture  distribution.  The  computed  pattern  obtained  by  applying  the 
synthesis  procedure  is  shown  in  Figure  7.  The  set  of  measured  data  points 
employed  (as  obtained  from  the  measured  pattern  of  Figure  6b)  are 
shown  on  the  pattern  for  comparison.  The  computed  values  of  the  12  line 
sources  from  which  the  synthesized  pattern  was  computed  are  also  shown, 
as  is  the  pattern  error  value  obtained  after  solution  convergence  was 
attained  (i.e.,  after  50  iterations,  for  the  data  shown  in  Figure  7). 

The  computed  pattern  of  Figure  7,  obtained  by  use  of  12  sources  to 
model  the  aperture,  agrees  very  closely  with  the  measured  pattern  from 
which  it  was  synthesized  (i.e..  Figure  6b).  In  order  to  demonstrate  how 
the  quality  of  such  solutions  may  degrade  if  too  few  sources  are  used  to 
model  the  aperture  distribution,  the  synthesis  procedure  was  applied  again 
to  the  same  problem,  but  with  only  7 sources  employed  to  model  the 
aperture  distribution.  The  resulting  solution  is  shown  in  Figure  8.  The 
relatively  poor  agreement  between  the  computed  and  measured  patterns  of 
Figure  8 and  Figure  6b,  respectively,  is  readily  apparent,  and  is  thought 
to  be  typical  of  the  solution  difficulties  which  may  arise  if  too  few 
sources  are  employed  to  model  the  aperture.  As  shown  by  the  previously 
considered  examples,  this  situation  can  be  avoided  by  employing  source 
spacings  on  the  order  of  a quarter  wavelength. 
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Figure  5.  Computed  solution  for  the  antenna  of  Figure  3,  obtained 
by  using  fewer  data  points  from  the  low  signal  level 
regions  of  the  measured  pattern.  It  may  be  noted  that 
the  solution  obtained  is  virtually  identical  to  the 
converged  solutions  of  Figures  4d  and  4e. 


- MEASURED 
. DATA  POINTS 


Computed  source  values  and  corresponding  pattern  for  the  antenna 
of  Figure  6.  As  shown  in  the  insert,  12  elemental  line  sources 
were  used  to  model  the  antenna  aperture. 


The  example  solutions  presented  thus  far  have  demonstrated  that  the 
synthesis  procedure  provides  a useful  method  for  computing  the  aperture 
distribution  of  an  antenna,  given  its  measured  radiation  pattern  and 
basic  physical  dimensions.  As  discussed  in  the  Introduction,  the  ultimate 
intended  purpose  of  the  procedure  is  to  be  able  to  employ  an  antenna's 
computed  aperture  distribution  to  predict  the  same  antenna's  radiation 
pattern  performance  in  a different  structural  environment  (e.g.,  when  the 
antenna  is  mounted  on  an  aircraft  fuselage).  As  a demonstration  of 
extending  the  analytical  solution  in  this  manner,  the  antenna  geometry 
of  Figure  3 was  modified  for  further  testing.  To  simulate  its  intro- 
duction into  a new  structural  environment,  an  adjustable  ground  plane  was 
added  to  the  original  antenna  geometry,  as  depicted  in  Figure  9.  Then, 


Figure  9.  Modified  version  of  the  antenna  geometry  of  Figure  3. 
Measured  and  computed  patterns  pertaining  to  this 
geometry  are  presented  in  Figure  10. 

using  the  computed  aperture  distribution  for  the  original  geometry  (as 
listed  in  Figure  4d)  together  with  a two-plate  diffraction  theory  model 
of  the  modified  antenna's  ground  plane  structure,  radiation  patterns  were 
computed  for  the  new  geometry,  and  compared  with  pattern  measurements  of 
the  new  geometry.  The  computed  and  measured  patterns  obtained  are 
compared  in  Figure  10,  for  several  selected  positions  of  the  adjustable 
ground  plane.  It  may  be  noted  that  although  the  computed  and  measured 


patterns  are  not  in  perfect  agreement,  the  extent  of  agreement  is  thought 
to  be  very  good,  especially  in  view  of  the  fact  that  the  computed 
solution  originated  from  a different  measured  pattern. 

A basic  question  of  validity  in  applying  the  above  technique  hinges 
on  the  validity  of  the  assumption  that  the  antenna's  aperture  distribution 
remains  essentially  invariant  when  the  antenna  is  introduced  into  a new 
structural  environment.  The  comparison  between  computed  and  measured 
patterns  presented  in  Figure  10  provides  strong  evidence  that  this 
assumption  is  valid  for  the  modified  geometry  considered.  As  a more 
direct  check  of  the  invariance  of  the  aperture  distribution,  the 
synthesis  procedure  was  applied  directly  to  both  the  original  antenna 
geometry  of  Figure  3 and  the  modified  geometry  of  Figure  lOd  (i.e.,  the 
case  of  a 5"  planar  extension  added  to  the  ground  plane  of  the  original 
antenna).  The  resulting  computed  aperture  distributions,  obtained  for 
each  geometry  after  50,  60,  70,  and  80  iterations  of  the  respective 
solutions,  are  listed  in  Table  I.  Both  solutions  have  essentially 
converged  after  40  iterations,  and  the  indicated  changes  in  either  of 
the  respective  solutions  after  50,  60,  70,  or  80  iterations  are  simply 
minor  oscillatory  adjustments  inherent  in  the  solution,  as  it  progresses. 
By  comparing  the  aperture  distributions  computed  for  the  respective 
antenna  geometries  after  a given  number  of  iterations  (e.g.,  after 
80  iterations  of  each  solution),  it  may  be  seen  that  the  distributions 
computed  for  the  respective  geometries  are  in  very  close  agreement.  In 
fact,  since  the  extent  of  agreement  between  solutions  is  roughly 
conmens urate  with  that  obtained  within  each  solution  as  the  numbe*  of 
solution  iterations  changes,  the  distributions  computed  for  the  cwo 
respective  geometries  could  be  viewed  as  being  virtually  identical. 


IV.  CONCLUSIONS 

A synthesis  procedure  is  presented  which  provides  an  efficient  means 
for  numerically  computing  the  aperture  distribution  of  an  antenna,  based 
on  its  specified  radiation  pattern  characteristics.  The  viewpoint 
adopted  here  differs  substantially  from  that  of  the  traditional 
synthesis  problem,  however,  in  that  the  specified  pattern  is  considered 
to  consist  of  measured  pattern  data  obtained  from  an  existing  well-defined 
antenna  geometry.  As  applied,  the  solution  of  importance  is  the  computed 
aperture  distribution  for  the  given  antenna.  Subsequently,  the  computed 
aperture  distribution  may  be  used  to  predict  (i.e.,  compute)  the  radiation 
pattern  performance  of  the  given  antenna  in  a different  structural 
environment  as,  for  example,  when  mounted  on  an  aircraft.  The  anticipated 
utility  of  such  a procedure  when  used  in  conjunction  with  available 
computer-aided  design  techniques,  provided  the  primary  motivation  for 
this  investigation. 

The  methods  presented  are  applicable  to  antennas  having  either 
continuous  or  discrete  apertures  mounted  in  a finite  or  infinite  ground 
plane.  Since  the  given  pattern  data  are  assumed  to  be  specified  in 
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RELATIVE  POWER  (dB)  RELATIVE  POWER  (dB) 


Cb)  120°  angle  between  plates 

Figure  10.  A specific  example  of  how  the  results  of  the  synthesis  procedure 
can  be  used  for  pattern  prediction.  The  computed  patterns  shown 
are  for  the  modified  antenna  of  Figure  9,  and  were  obtained  by 
using  the  computed  source  distribution  for  the  unmodified  antenna 
of  Figure  3 (i.e.,  the  source  values  of  Figure  4d).  A two-plate 
model  is  used  to  simulate  the  structure.  The  measured  patterns 
shown  were  obtained  directly  from  the  antenna  geometry  depicted 
in  Figure  9. 
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ANGLE  6 ( DEGREES  ) 


(c)  150°  angle  between  plates 
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ANGLE  6 (DEGREES) 

(d)  180°  angle  between  plates 


Figure  10  (Continued). 


TABLE  I 

Source  Values  Obtained  After  50,  60,  70  and  80  Iterations 
For  The  Antennas  of  Figure  3,  and  of  Figure  9 With  a=180° 
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magnitude  only,  a numerical  iteration  process  is  required  to  evolve  the 
unspecified  phase  information.  In  cases  where  both  pattern  magnitude  and 
phase  data  are  available,  the  solution  for  the  desired  aperture  distri- 
bution may  be  obtained  directly  without  using  the  iteration  process. 

In  addition  to  principal  plane  (e.g.,  E-  and  H-plane)  pattern  data, 
application  of  the  synthesis  procedure  requires  that  the  physical  dimensions 
of  the  antenna  aperture  and  the  dimensions  of  the  ground  plane  in  which 
it  is  mounted  (if  any),  be  given.  For  antennas  mounted  in  finite  ground 
planes,  the  Geometrical  Theory  of  Diffraction  (GTD)  is  used  to  account 
for  the  effects  of  the  ground  plane  edges.  The  antennas  specifically 
considered  in  this  investigation  were  composed  of  rectangular  apertures, 
flush  mounted  in  planar  rectangular  ground  planes.  For  the  apertures 
and  pattern  planes  considered,  the  aperture  distributions  could  be 
treated  as  one-dimensional.  The  basic  methods  presented  are  more 
generally  applicable,  however,  and  could  be  extended  to  treat  non-planar 
apertures  of  arbitrary  shape,  which  are  flush  mounted  in  curved  surfaces. 

During  the  course  of  the  investigation  the  synthesis  procedure  was 
applied  to  several  specific  antenna  geometries.  The  solutions  obtained, 
as  presented  in  Section  III,  serve  to  demonstrate  the  capabilities  of 
the  procedure.  They  also  serve  to  point  out  some  specific  techniques 
which  were  found  to  be  helpful  in  successfully  implementing  the  procedure. 

It  was  found,  for  example,  that  the  discrete  sources  employed  to  model 
the  aperture  distribution  should  be  spaced  about  a quarter  wavelength 
apart.  As  demonstrated  by  one  of  the  examples,  use  of  larger  source 
spacings  (e.g.,  spacings  of  A/2)  may  result  in  relatively  large  errors  in 
the  synthesized  pattern  (indicating  use  of  an  inadequate  aperture 
representation).  Additional  helpful  techniques,  relating  to  the  most 
appropriate  selection  of  data  points  from  a given  antenna  pattern,  were 
also  deduced  from  the  antenna  geometries  considered.  When  the  antenna 
aperture  is  mounted  in  a finite  ground  plane,  it  was  found  that  the  most 
efficient  selection  of  data  points  is  obtained  when  all  of  the  points 
are  chosen  from  the  illuminated  half  of  the  radiation  sphere  (i.e., 
from  the  aperture  side  of  the  ground  plane).  The  reason  for  this  is  that 
the  diffracted  fields  on  the  shadow-side  of  the  ground  plane  (i.e.,  the 
total  fields  in  this  region)  are  only  indirectly  related  to  the  aperture 
fields  (through  the  field  incident  on  the  edge)  and  these  fields  are 
adequately  represented  by  their  inclusion  in  the  illuminated  half-space 
fields. 

Aside  from  the  data  selection  consideration  just  noted,  one  has  a 
considerable  degree  of  latitude  in  selecting  both  the  number  and  relative 
location  of  pattern  data  points.  When,  as  assumed  here,  the  pattern 
data  points  are  selected  from  a measured  pattern,  it  may  be  helpful  to 
select  fewer  data  points  from  the  low  signal  level  regions  of  the  pattern. 
This  technique,  which  was  demonstrated  in  the  previous  section,  simply 
acknowledges  the  fact  that  errors  in  measured  signal  levels  tend  to  be 
larger  in  the  low  level  regions  of  the  pattern.  Choice  of  fewer  inaccurate 
data  values  often  yields  as  accurate  a solution  as  can  be  obtained,  with 
an  expenditure  of  considerably  less  computation  time.  Finally,  it  has 
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been  observed  that  solution  difficulties  may  result  if  the  true  angular 
orientation  of  the  measured  pattern  data  is  not  accurately  specified. 
Specifically,  if  the  angular  orientation  of  the  pattern  data  is 
incorrectly  specified  by  even  a small  amount  (e.g.,  1 or  2 degrees), 
solutions  were  found  to  converge  rapidly,  but  the  resulting  synthesized 
pattern  typically  exhibited  a relatively  large  pattern  error.  It  was 
found  that  this  difficulty  could  be  overcome,  however,  by  rotating  the 
entire  data  set  by  a small  amount,  and  repeating  the  solution  until  the 
orientation  yielding  the  smallest  pattern  error  is  determined.  This 
can  only  be  done  for  small  angular  error  in  the  alignment  of  the  ground 
plane.  Consequently,  one  should  attempt  to  obtain  the  best  possible  set 
of  measurements  for  the  antenna  under  consideration.  It  is  obvious  that 
one  can  not  start  with  patterns  with  large  measurement  errors  and  hope  to 
synthesize  an  error-free  aperture  distribution. 
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APPENDIX  A 


This  appendix  contains  a brief,  user  oriented  description  and 
statement  listing  of  the  computer  program  developed  as  a part  of  the 
investigation.  The  program  presented  was  used  to  obtain  the  computed 
results  presented  in  Section  III  of  the  text.  As  presented,  the  program 
is  written  for  use  on  an  interactive  system  permitting  teletype  input 
and  output.  Provision  is  also  included  for  visual  cathode  ray  tube  (CRT) 
display  of  the  computed  patterns,  generated  at  specified  intervals  as 
the  iterative  solution  progresses.  Programming  statements  associated  with 
these  capabilities  and  other  minor  provisions  which  may  not  be  con- 
veniently implemented  on  other  facilities,  are  noted  with  asterisks  in 
the  Drogram  listing.  Comment  statements  (i.e.,  statements  preceded  by 
C!!!)  are  distributed  throughout  the  program  to  describe  the  analytical 
function  of  various  statement  sequences. 

As  an  aid  in  discussing  the  program,  its  major  sections  are  depicted 
in  block  form  in  Table  II.  A brief  description  of  the  major  functions 
performed  within  each  of  the  letter-designated  blocks  is  given  as  follows: 

A - This  block  of  statements  (lines  21-52)  inputs  the  antenna  geometry, 
the  integer  number  of  sources  (NS)  chosen  to  model  the  aperture,  the 
integer  number  (ITMX)  of  iterations  desired,  and  the  list  of  data  values 
selected  from  the  antenna's  measured  pattern.  The  antenna  geometry  is 
specified  by  inputting  the  dimensions  SMIN,  SMAX,  and  W,  expressed  in 
wavelengths.  Definition  of  these  quantities  is  depicted  in  Figure  11. 


Figure  11.  Geometry  defining  various  antenna  and  pattern 
parameters  used  in  the  computer  code. 
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TABLE  II 


START 


END 
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The  pattern  data  values  are  defined  by  integer  number  (M),  angle  PH(M),  and 
magnitude  HD(M),  as  depicted  in  Figure  11.  This  data  is  read  in  free- 
format,  from  a separate  file.  The  specific  list  of  pattern  data  values 
used  in  the  example  of  Figure  4 (Section  III)  is  shown  in  Table  III. 

Note  that  the  first  number  to  be  read  is  the  total  integer  number  of  data 
values  contained  in  the  list. 

B - This  group  of  statements  (lines  53-58)  computes  the  individual 
elements  of  the  matrix  [TmnL  by  making  repeated  calls  to  the  function 
subroutine  HPLN.  This  subroutine  computes  the  total  field  (including 
diffraction  from  ground  plane  edges)  of  each  elemental  source  at  each 
field  point  location,  and  is  listed  following  the  main  program. 

C - This  sequence  of  statements  (lines  72-98)  sets  up  a CRT  display, 
consisting  of  a polar  grid  and  plot  of  the  selected  measured  data 
points.  Later,  in  the  iteration  loop,  a plot  of  the  computed  pattern 
is  added  to  the  display  for  the  purpose  of  visual  inspection  of 
solution  progress.  The  pattern  is  recomputed  for  display  purposes  at 
intervals  specified  by  the  integer  variable  NPLTMX  (read  in  statement  98). 

A value  of  5 for  NPLTMX,  for  example,  would  result  in  a visual  display 
update  of  the  computed  pattern  after  every  5th  iteration.  (If  desired, 
this  entire  section  may  be  bypassed  or  deleted,  except  for  specification 
of  a value  for  NPLTMX,  which  is  used  elsewhere  in  the  program  to  write 
updated  lists  of  the  computed  source  values.) 

D - This  group  of  statements  (lines  109-213)  iteratively  recomputes  the 
solution,  until  a convergent  result  is  obtained.  The  source  and  pattern 
error  values  computed  after  each  iteration  are  written  (the  error  values 
are  written  as  teletype  output)  at  intervals  specified  by  NPLTMX.  The 
computed  pattern  is  also  generated  and  displayed  on  the  CRT  at  the  same 
specified  interval,  as  discussed  above.  If  a visual  display  cannot  be 
used,  statements  160-172  should  be  bypassed  or  deleted. 

E - This  block  of  statements  (lines  217-226)  recomputes  the  total 
pattern  after  completion  or  termination  of  the  iteration  loop.  The 
computed  pattern  values  are  stored  in  a binary  file  (named  SYNSIS)  for 
later  use. 

F - This  group  of  statements  (lines  229-263)  initiates  an  optional 
polar  pen  plot  of  the  final  computed  pattern.  The  measured  pattern  data 
points  used  as  input  in  the .procedure  are  also  plotted  on  the  same  grid, 
for  comparison  with  the  computed  pattern.  The  plots  presented  in 
Section  III  (e.g.,  see  Figure  4)  were  obtained  by  using  this  plot  section. 
This  group  of  statements  may  be  bypassed  or  deleted,  if  desired. 

The  above  summary  completes  the  discussion  of  the  main  synthesis 
program,  which  extends  from  statements  1-265  in  the  accompanying  program 
listing.  Several  subroutines  used  in  the  main  program  are  also  included 
in  the  listing  shown.  Subroutines  GRID  and  POLPLT  are  used  in  the  various 
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plot  sections  of  the  main  program.  The  complex  function  subroutine  HPLN 
is  used  to  compute  the  source  fields  in  the  main  program.  Finally, 
subroutines  DI  and  FRNELS  are  diffraction  coefficient  and  fresnel  integral 
subprograms  required  by  the  source  subroutine  HPLN.  Two  standard 
routines  used  by  the  program  are  not  listed.  One  of  these  is  the  sub- 
routine PLOT,  which  is  a standard  CALCOMP  plot  routine.  The  other 
subprogram  not  specifically  listed  is  subroutine  CGELG.  This  is  a 
standard  IBM  program  for  solving  simultaneous  equations,  which  was 
modified  slightly  to  treat  complex  arithmetic. 
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TABLE  III 

LIST  OF  MEASURED  PATTERN  DATA  VALUES  USED  AS  COMPUTER  CODE 
INPUT  TO  OBTAIN  THE  COMPUTED  SOLUTIONS  OF  FIGURE  4 


91 

0.,- 

2. 

4.  ,- 
6. 

8.  ,- 
10.  , 
12.  , 
14.  , 
16.  , 
18.  , 
20.  , 
22.  , 
24.  , 
26.  , 
28.  , 
30.  , 
32.  , 
34.  , 
36.  , 
38.  , 
40.  , 
42.  , 
44.  , 
46.  , 
48.  , 
50.  , 
52.  , 
54.  , 
56.  , 
58.  , 
60.  , 
62.  , 
64.  , 
66.  , 
68.  , 
70.  , 
72.  , 
74.  , 
76.  , 
78.  , 
80.  , 
82.  , 
84.  , 
86.  , 
88.  , 
90.  , 


■13.3 
■12.25 
■11.25 
10.15 
•9.15 
-8.1 
-7.3 
-6.45 
-5.8  5 
-5.3 
-4.95 
-4.7  5 
-4.7 
-4.75 
-4.9 
-5.15 
-5.25 
-5. 
-4.45 
-3.8 
-3.15 
-2.75 
-2.6  5 
-2.7 
-2.9 
-3.0  5 
-3. 
-2.6 
-2.25 
-2.1 
-2.25 

- 2.6  5 
-2.9 

- 3. 

- 2.8  8 

- 2.8  5 

- 3. 
-3.7  5 
-4.35 
-4.8 
-4.8  5 
-5.15 
-5.8 
-6.6 
-7.3 
-7.7 


92.  , 

-7 

.85 

94.  , 

-8 

. 3 

96.  , 

-9 

• 

98.  , 

-9 

.8 

100. 

10 

• 

102. 

10 

• 

104. 

10 

. 3 

106. 

11 

. 2 

108. 

12 

. 

110. 

12 

. 3 

112. 

t — 

12 

. 3 

114. 

11 

.8 

116. 

12 

.0  5 

118. 

13 

. 2 

120. 

14 

.3 

122. 

, - 

14 

.8  5 

124. 

14 

.9 

126. 

13 

.75 

128. 

, - 

13 

. 15 

130. 

, - 

13 

.25 

132. 

/ — 

13 

.8  5 

134. 

/ — 

14 

.8 

136. 

,- 

16 

. 

138. 

, - 

16 

.3 

140. 

/ — 

16 

.28 

142. 

/ — 

15 

.6  5 

144. 

t — 

15 

• 

146. 

, - 

15 

• 

148. 

15 

.15 

150. 

/ — 

14 

.6 

152. 

14 

.6 

154. 

15 

.5 

156. 

16 

• 

158. 

, - 

15 

.8  5 

160. 

15 

.8 

162. 

, - 

16 

.4 

164. 

16 

.6 

166. 

16 

.8 

168. 

17 

. 3 

170. 

18 

.05 

172. 

18 

.7 

174. 

18 

. 8 

176. 

, - 

19 

.15 

178. 

19 

.85 

180. 

20 

.5 

1 OPTIONS  DP 

2 OPTIONS  32K 

3 INCLUDE  CGELGB  , 3 39  OM  ;FLSBS  , 4 0 4C 

4 COMPLEX  HS (120,4 0) , TS, HITM( 120) , F(40 ) , T(40, 4 0) 

5 2,  C(  16  00  ),  HIT (361)  , HPLN,  CJ,  FF  (40)  , HLAST  (120  ) 

6 COMPLEX  TOP, FAC 

7 DIMENSION  IFI  LE  ( 2 ) , I US  ER  ( 2)  , PH  ( 12  0)  , H D ( 12  0)  , 

8 2PHAR(  120)  , IBUF  (5  12) 

9 LOGICAL  LOG , LV  IS  D,  LDB  , LSO  R,  L PH  AS  E 

10  COMMON/CONST/TOP 

11  COMMON/S  TO RE/T,C,HS 

12  COMMON/ P LAT/W , LSO R 

13  EXTERNAL  CRTP LT , P LOT 

14  LSOR=.  FALSE. 

15  PI=3. 14159 265 

16  TPI=2.*PI 

17  DPR=1 8 0 . / P I 

18  RPD=1./DPR 

19  CJ=(0.,1.) 

20  TOP=-CEXP(-CJ*0.25*PI) 

21  WRITE (8,1) 

22  1 FORMAT  ( ' MIN  AND  MAX  DISTANCE  OF  APERTURE  FROM  ED  GE=  ? ') 

23  READ  (8  , - ) SMIN,  SMAX 

24  WRITE (8,6  ) 

25  6 FORMAT  ( * PLATE  WIDTH  IN  WAVELE  NGTH  S=?  ') 

2 6 READ  (8,-)  W 

27  WRITE  (8,2) 

28  2 FORMAT  ( ' NUMBER  OF  SOURCES  TO  SIMULATE  APERTURE=  ? ') 

29  READ  (8,-)  NS 

*30  CALL  ESC ($9898) 

31  WRITE  (8,30) 

32  30  FORMAT  ( ' NUMBER  OF  ITERATIONS  DE  SI  RED=?  ') 

33  READ  ( 8,  - ) ITMX 

34  Cl  11  COMPUTE  SOURCE  SPACING. 

35  DELS=  (SMAX -SMIN)  / NS 

36  WRITE  (8,23) 

37  23  FORMAT  ('  IS  THE  PATTERN  DATA  TO  BE  INPUT  IN  DECIBELS 

38  2 (T  OR  F)  ?'  ) 

39  READ  (8,-)  LDB 

40  WRITE (8,3) 

41  3 FORMAT  ( * INPUT  FILE  NAME  OF  PATTERN  DATA  POINTS:  ' ) 

*42  CALL  RDFLNM(  IFILE  , I US  ER) 

*43  CALL  ASSIGN  (IFILE  , IUSER,7  ) 

44  READ (7,-)  MPMX 

45  READ  ( 7 , - ) ( PH  (M)  , HD  (M)  , M=  1,  MPMX  ) 

46  Clll  CONVERT  DECIBEL  PATTERN  VALUES  TO  FIELD  VALUES. 

4 7 IF  (.NOT.  LDB)  GO  TO  24 

48  DO  25  M=1 , MPMX 

49  25  HD  (M)  =1  0 . * * ( . 0 5*H  D (M)  ) 
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50 

24 

51 

52 

5 

53 

C!  ! i 

54 

55 

56 

57 

58 

4 

59 

C 1 l l 

60 

C ! ! ! 

61 

62 

63 

64 

10 

65 

C 1 t 1 

66 

67 

68 

69 

70 

12 

71 

11 

*72 

*73 

22 

*74 

*75 

*76 

CM! 

*77 

*78 

*79 

*80 

*81 

*82 

*83 

20 

*84 

*85 

*86 

*87 

*88 

*89 

*90 

*91 

*92 

21 

*93 

CM! 

*94 

C ! ! ! 

*95 

*96 

44 

*9? 

38 

CONTINUE 
DO  5 M=1,MPMX 
PH ( M) =180.- PH (M ) 

COMPUTE  MATRIX  ELEMENTS  FOR  SOURCE  MATRIX  T. 

DO  4 N=1,NS 
S=SMIN+  (N-  0.5  ) * D E LS 
DO  4 M= 1 , MPMX 
PHR=PH (M) * RPD 
HS(M,N)=HPLN(S,PHR) 

SET  INITIAL  ZERO  PHASE  AND  COMPUTE  NORMALIZATION  FOR 
ERROR  CRITERION  (ER) 

HOS=0 . 

DO  10  M=1 / MPMX 
PHAR(M) =0. 

HOS=HOS+HD (M) *HD(M) 

COMPOTE  MATRIX  ELEMENTS 
DO  11  NI=1,NS 
DO  11  NJ= 1, NS 
TS=  (0.,0. ) 

DO  12  M=1 » MPMX 

TS=TS+C  ON  JG  ( HS  (M,  NI)  ) *HS(M,NJ) 

T (NI,  NJ)  =T  S 
WRITE  (8,2  2) 

FORMAT  ( ' IS  CRT  DISPLAY  TO  BE  USED(T  OR  F)?') 

READ ( 8 , - ) LVISD 

IF  ( .NOT.  LVISD)  GO  TO  44 

SET  UP  CRT/PLOT  GRID, AND  DESIRED  VALUES 
CALL  CRTON 

CALL  CRTPTS ( IB  UF , 1 0 , 3 ) 

CALL  BRPLOT  (C  RT  PLT) 

CALL  GRID 
HMAX=0. 

DO  20  M= 1 , MPMX 

IF(HD(M) .GT.HMAX)  HMAX=HD (M) 

DO  21  M= 1 , M PMX 

HPLT=20 . *A  LOG 1 0 ( HD (M) /HMAX ) 

IF(HPLT  .LT.-  40.  ) HPLT=-40. 

HMX= (40.+HPLT) *2.5/40. 

PHR=PH (M) *R  PD 
PHR=PHR- 0. 5 *P I 
XPLT=HMX*S  IN  (PH  R) 

YPLT=HMX*COS  ( PH  R) 

CALL  SYMBOL  ( XPLT  ,Y  PLT,  0 . 0 25, 1 1,  0 . ,-  1) 

GET  ADDRESS  OF  GRID  AND  DESIRED  VALUE  TO  ERASE  PLOT 
DATA. 

CALL  INFO  (XPLT,  Y PLT  , X O,  YO,  FACT  , I TE  N,  I AD  R) 

WRITE  (8, 3 8) 

FORMAT  ( * INCREMENT  DATA  PLOTTED  BY  = ? ') 
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98 

READ  (8,-)  NPLTMX 

99 

NPLT=NPLTMX 

100 

MI=1 

101 

MF=MPMX 

10  2 

MS=  1 

103 

TESTN=1.  j 

104 

DO  821  M=MI,MF,MS 

105 

821 

BLAST (M) = (0.  ,0. ) 

106 

LPHASE=  . TRUE  . 

107 

DELT= . 0 4 

103 

DELTN=-DELT 

109 

C ! ! 1 

ITERATE  ON  RESULTS 

110 

DO  31  IT= 1 , IT  MX 

111 

NPLT=NPLT+  1 

112 

C J * i 

COMPUTE  COMPLEX  FIELD  POINT  VALUES. 

113 

DO  32  M=MI , MF , MS 

114 

HITM  ( M) =HD (M) *C  EXP (C J *P  HAR(M)  ) 

115 

C ! ! ! 

COMPUTE  UPDATED  FIELD  POINT  VALUES 

116 

C ! i ! 

BASED  ON  INCREMENTAL  PHASE  CHANGE  j 

117 

Ciii 

BETWEEN  LAST  AND  PRESENT  ITERATION. 

118 

HITM(M)  =H  LA  ST  (M  ) +TESTN*  ( HIT  M(M)  -HLAST(M)  ) 

119 

32 

HLAST  (M)=HITM(M) 

120 

Ciii 

COMPUTE  SOURCE  VECTOR. 

121 

DO  33  N=1,NS 

122 

% 

z 

II 

o 

• 

O 

• 

123 

DO  33  M=MI,MF,MS 

124 

3 3 

F ( N ) =F  (N  ) + H IT  M(M ) *CONJ G ( HS  (M,  N)  ) 

125 

DO  34  N J= 1 / NS 

126 

DO  34  NI= 1 / NS 

127 

NN=  (NJ-  1)  * N S+  N I 

128 

34 

C (NN)  =T  (NI,  NJ) 

129 

C I ! • 

SOLVE  SIMULTANEOUS  EQUATIONS. 

130 

CALL  CGELG(F,C,NS,l,l.E-6,  IER) 

131 

IF(IER.NE.O)  WR ITE (8,35)  IER 

132 

35 

FORMAT  ( * ERROR  IN  GELG  : IE  R=' , I 2) 

133 

IF(LPHASE)  GO  TO  521 

134 

DO  534  NJ=1,NS 

135 

DO  534  NI=1 / NS 

136 

NN=  (NJ-  1)  * N S+  N I 

137 

5 34 

C ( NN)  =T  (N  I,  NJ) 

138 

cut 

COMPUTE  UPDATED  FIELD  POINT  DIFFERENCE  VALUES 

139 

cut 

BASED  ON  INCREMENTAL  MAGNITUDE  CHANGE  BETWEEN 

140 

C!  U 

ITERATION  JUST  COMPLETED  AND  PREVIOUS  ITERATION. 

141 

DO  532  M=MI,  MF  , MS 

142 

TS=  (0  . ,0.  ) 

143 

DO  55  2 N=  1 , NS 

144 

55  2 

TS=TS+F  (N) *H  S (M ,N ) 

145 

PHARN=CATAN2  (TS) 

146 

HMN=CABS  (TS) 

35 


r* 

1 4 7 

532 

HITM(M)  =HD(M)  *C  EX  P(  C J*  P HARN  ) -HMN*C  EX  P (C  J*P  HAR  (M)  ) 

148 

C ! i i 

COMPUTE  INCREMENTAL  CHANGE  IN  SOURCES 

149 

C ! i i 

DUE  TO  CHANGE  IN  FIELD  POINT  VALUES. 

150 

DO  533  N=l, NS 

151 

FF(N)= (0. ,0. ) 

152 

DO  533  M=MI,MF,MS 

153 

533 

FF(N)=FF(N)  +H  ITM(M)  *CONJG(HS  (M,N)) 

154 

CALL  CG EL G(FF,C,NS,l,l.E-6,  IER) 

155 

IF(IER.NE.O)  WRITE  (8,35)  IER 

156 

c • i • 

ADD  INCREMENTAL  CHANGE  IN  SOURCES  TO 

157 

C ! ! 1 

PREVIOUS  SOURCE  VALUES. 

158 

DO  522  N=l, NS 

159 

522 

F(N)=F  (N)+FF(N) 

* 1 60 

521 

IF ( .NOT.LVISD)  GO  TO  41 

j 

* 161 

C ! ! I 

COMPUTE  TOTAL  PATTERN. 

* 162 

IF  (NPLT  .LT.  NPLTMX)  GO  TO  41 

* 163 

DO  40  I P=  0 , 3 60 

* 164 

I D= I P+ 1 

* 165 

PH R=  I P*R  PD 

* 166 

HIT ( ID ) = (0.  ,0.  ) 

* 167 

DO  40  N= 1 / N S 

* 168 

S=SMIN+  (N-0.5  ) *D  ELS 

* 169 

40 

HIT  (ID)  =HIT  (ID)  +F  (N)  * H PLN  (S  ,PHR) 

* 170 

CALL  ERASE  (IADR) 

* 171 

CALL  CRTPLT  (4 . 2 5,  4 . 25,-  3) 

* 172 

CALL  POLPLT(HIT) 

173 

41 

FNS  = 0 . 

174 

HMS=0 . 

175 

HHS  = 0 . 

176 

DO  50  N= 1 , N S 

177 

50 

FNS=F  NS+C  AB3  ( F (N  ) ) * *2 

178 

DO  51  M=1 , MPMX 

179 

TS=  (0.  ,0.  ) 

180 

DO  52  N= 1 , NS 

181 

52 

T S=T  S+F  (N)  * H S (M  , N ) 

182 

PHAR(M) = C AT  AN  2 (TS) 

183 

HM=CAB3 ( TS ) 

184 

HMS  = HMS+HM*HM 

185 

51 

HHS  = HHS+  (HM-HD(M)  ) * *2 

186 

ER=HHS/HOS 

187 

QR=MPMX  *F  NS/HMS 

188 

WRITE  ( 8,  - ) I T,  ER  , QR 

189 

IF(IT.LE.5  ) DE  LT  N=  -D  ELT 

190 

DELTN=DE  LT  N+DELT 

191 

IF  ( ER  . GT.  ERLAST)DELTN=DELTN*0.5 

192 

E RLAS  T=E R 

193 

IF  ( IT.  GE  . 5 ) L PH  AS  E=  . FALSE  . 

194 

TESTN=1 .+DELTN 

195 

IF(TESTN.  GT.  2 . ) TE STN=  2 . 

196 

IF  (NPLT  .LT.  NPLTMX)  GO  TO  10  31 

197 

WRITE (8,7 23)TE^TN 

198 

723 

FORMAT  ( ' WEIGHTING  FACTOR=  ',F10.5) 

36 

199 

200 

201 

202 

203 

204 

205 

206  643 

207 

208 

209 

210  46 

211  1031 

212 

213  31 
*214  9898 

215  C ! ! 1 

216  C i ! ! 

217 

218 

219 

220 
221 
222 

223  45 

224  C ! ! ! 

225  C ! ! i 
*226 
*227 
*228 

*229  C ! 1 ! 
*230  C ! ! ! 
*231  9899 
*232 
*233 
*234 
*235 
*236 

*237  9376 

*238 

*239 

*240 

*241 

*242  720 

*243 

*244 

*245 

*246 

*247 


WRITE ( 6 , - ) IT 
ICMAX=0 
CURMAX=  0. 

DO  643  1=1, NS 

IF (CABS ( F( I) ) . LT .CURMAX)  GO  TO  643 
CURMAX=CABS(F (I ) ) 

ICMAX=I 
CONTINUE 
DO  46  1=1, NS 

CUR M=  C A BS  (F  ( I)  /F  ( ICMAX)  ) 

CURP=DP  R*C  AT AN 2 (F  (I)/F  (ICMAX) ) 

WRITE  (6,-)  I , C URM,  CU  RP 

IF  (NPLT  .GE.  NPLTMX  ) NPLT=0 

IF  ( IT  . EQ.  1 ) N PLT=  1 

CONTINUE 

CALL  ESC(  $99  99  ) 

COMPUTE  TOTAL  PATTERN  AFTER  COMPLETION 
OR  TERMINATION  OF  ITERATION  LOOP. 

DO  45  I P=  0 , 3 60 
ID=IP+1 
PHR=  I P*  R PD 
HIT ( ID ) = (0.  ,0.  ) 

DO  45  N=1,NS 
S=SMIN+  (N- . 5 ) *D  ELS 

HIT (ID)=HIT  (ID) +F (N) *H  PLN (S ,PHR) 

WRITE  COMPUTED  PATTERN  VALUES  INTO 
PERMANENT  FILE  NAMED  SYNSIS. 

CALL  ASSIGN (6HSYNS IS, 0. ,1) 

WRITE  ( 1 ) HIT 
WRITE  ( 8,  9 899) 

PLOT  COMPUTED  PATTERN  AND  MEASURED  PATTERN 
VALUES  FROM  WHICH  PATTERN  WAS  COMPUTED. 
FORMAT  ('  PEN  PLOT  DE  SI  RE  D ( T OR  F)?') 

READ  (8,-)  LOG 
IF ( . NOT .LOG ) GO  TO  9999 
IF  ( .NOT  .LVISD)  GO  TO  93  76 
CALL  PLOT  (0  . ,0  . ,-  3) 

CALL  PLOT ( 0 . ,0.  ,999  ) 

CALL  PLOTS(  IBUF  ,5  12, 3 ) 

CALL  BRPLOT ( PLOT ) 

CALL  GRID 
HMAX=0 . 

DO  7 20  M= 1 , MP MX 

IF (HD (M)  .GT .HMAX)  HMAX=H  D (M ) 

DO  721  M=1 , MPMX 

HPLT=  2 0 . *A  LOG1 0 (HD (M) /HMAX) 

IF(HPLT .LT.- 40. ) HPLT=-40. 

HMX=  (40.+  HPLT)*2.5/40. 

PHR=PH (M) * R PD 


g— ■ BH  " 


*248 

PHR=PHR- 0.5  *PI 

*249 

XPLT=HMX*S  IN(PHR) 

*250 

YPLT=HMX*COS(  PH  R) 

*251 

721 

CALL  SYMBOL(XPLT,YPLT,  0.025,11,0.,-  1) 

*252 

CALL  POLPLT  (H  IT  ) 

*253 

CALL  PLOTS ( IBUF ,512,3 ) 

*254 

LSOR= . TRUE . 

*255 

DO  786  IP=  0 , 3 6 0 

*256 

I D=  I P+ 1 

*257 

PHR=IP*RPD 

*258 

HIT ( ID ) = ( 0 . ,0. ) 

*259 

DO  786  N=l, NS 

*260 

S=SMIN+ (N-.5  ) * D EL  S 

*261 

786 

HIT(ID)=HIT(  ID  ) +F  (N)  *HPLN(S  ,PHR) 

*262 

CALL  GRID 

*263 

CALL  POL 9LT ( HIT ) 

264 

9999 

CALL  EXIT 

265 

END 

266 

C!  ! ! 

****  END  OF  MAIN  PROGRAM  **** 

267 

SUBROUTINE  GRID 

268 

DATA  PI, TPI ,D PR/3 . 14159  26 5,  6 . 2 8 318  53, 5 7. 2 95  77  95  8/ 

269 

RP=  2.5 

270 

CALL  PLOT  (4.25,4.25,-  3) 

271 

C J J 1 

POLAR  GRID 

272 

DO  110  1=1,4 

273 

RG=RP*I/4 . 

274 

CALL  PLOT ( RG,  0 . ,3  ) 

275 

DO  110  J=0 ,360,2 

276 

ANG= J/D PR 

277 

XX=RG*C  OS ( ANG ) 

278 

YY=RG*S IN (ANG) 

279 

110 

CALL  PLOT(  XX,  Y Y,  2 ) 

280 

DO  112  1=1,6 

281 

ANG=  (I-  1)  *PI/6. 

282 

ANGS=ANG+PI 

283 

ANGF= ANG 

284 

IF(I.EQ.2* (1/2))  GO  TO  111 

285 

ANGS  = AN  G 

286 

ANGF=ANG+P I 

287 

111 

CO  NT  INUE 

288 

XX=RP*COS  ( ANG  S) 

289 

YY=RP*S  IN  (ANGS) 

290 

CALL  PLOT  (XX,  YY,  3 ) 

291 

XX=RP*COS  ( ANGF) 

29  2 

YY=RP*S  IN  (ANGF) 

293 

112 

CALL  PLOT(  XX,  YY,  2 ) 

294 

RETURN 

295 

END 

296 

SUBROUTINE  POLPLT(ET) 

297 

COMPLEX  ET (361) 

38 


298 

DATA  PI,  TPI, DPR/3. 141  59  26 5,  6 .2831853,57.29577958/ 

299 

RP= 2 . 5 

300 

EMX=0  . 

301 

DO  101  IP=0 ,360 

302 

I=IP+1 

303 

EM=CABS  ( ET  (I)  ) 

304 

IF  (EM.GT.EMX  ) EMX  =E  M 

305 

101 

CONTINUE 

306 

CM! 

PATTERN  PLOT 

307 

DO  120  I P=  0 , 3 6 0 

308 

I=IP+1 

309 

ETM=CABS  (ET  ( I)  ) /E  MX 

310 

IF  (ETM  .LT.  0 . 0 1)  ETM=0.01 

311 

RD=  2 0. *A  LOG1 0 (ETM) 

312 

IF(RD.LT.- 40.  ) RD=-40. 

313 

RD=RP* ( RD+40. )/40. 

314 

125 

CONTINUE 

315 

ANG= IP/D  PR 

316 

ANG=ANG-0 . 5 *PI 

317 

XX=RD*S  IN  (ANG) 

318 

YY=  RD*C  OS  (ANG) 

319 

I PEN=  2 

320 

IF(I.EQ.l)  IP  EN=  3 

321 

120 

CALL  PLOT  (XX,  YY,  IPEN) 

322 

CALL  PLOT (4 . 25,- 5.5 ,- 3) 

323 

13  0 

CONTINUE 

324 

CALL  PLOT (0.  ,0.  ,9  99  ) 

325 

RETURN 

326 

END 

327 

COMPLEX  FUNCTION  HPLN(S,PHR) 

328 

COMPLEX  C J , D C, H I 

329 

COMMON/P LAT/W ,LSOR 

330 

LOGICAL  LSOR 

331 

DATA  CJ/(0.  ,1.  )/ 

332 

DATA  PI,  TPI,  DPR/  3. 14159  265,6.2831853,57.29577958/ 

333 

HI=CEXP(-CJ*TPI*S  ) /SQRT(S) 

334 

HPLN=CEXP(CJ*T  PI*S*COS (PHR) ) 

335 

IF  (LSOR)  RETURN 

336 

CALL  DI (DC,  S,  PHR,1.  ,2.  ) 

337 

IF (PHR.GT. PI ) HP  LN=  ( 0 . , 0 . ) 

338 

HPLN=HPLN+HI*DC 

339 

H I=CEXP  (-CJ*TPI*  ( W-S  ) ) /SQRT(W-S) 

340 

PHR2  = PI-PHR 

341 

I F ( P HR  2 . LT . 0 . ) PHR2=TPI+  PHR2 

342 

CALL  DI (DC, W-S ,PHR2, 1 . , 2 . ) 

34  3 

HPLN=HPLN+HI*DC*CEXP(CJ*TPI*W*COS  (PHR)  ) 

344 

RETURN 

345 

END 

346 

SUBROUTINE  DI  (D  IR  , R , ANG  R,  SBO  , FN  ) 

39 


347  C ***  INCIDENT  (BET-PH-PHP)  OR  REFLECTED  (BET*PH+PHP)  **' 

348  C ***  PART  OF  WEDGE  DIFFRACTION  COEFFICIENT  *** 

349  COMPLEX  COM,  EX,  UPPI,  UNPI,FA,DIR 

350  COMPLEX  TOP 

351  COMMON/CONST/TOP 

352  DATA  PI , T PI , D PR/ 3 . 1 41  59  26  5, 6 . 2 8 31 8 53 , 5 7.  2 95  77  95  8/ 

353  DEM=  2.  *T  PI*F  N*S  BO 

3 54  COM=TOP/D  EM 

355  SQR=SQRT(TPI*R) 

356  DNS=  ( PI+ ANG  R)/(2.0*FN*PI) 

357  SGN=SIGN  (1.  ,DNS) 

358  N=IFIX  (ABS  (DNS) +0.5  ) 

359  DN=SGN*N 

36  0 A=ABS  (1.0+COS(  ANGR-  2.0*FN*P  I*D  N)  ) 

361  BOTL  = 2 . 0*SQRT  (ABS  (R*A)  ) 

362  EX=CEX  P( CMP  LX (0.0, TP I*R*A ) ) 

363  CALL  FRNELS  (C,S,BOTL) 

364  C=SQRT( PI/2.0 ) * ( 0 . 5 -C ) 

36  5 S=  SQRT(  P 1/2.0  ) * ( S-  0 . 5 ) 

36  6 F A=CMPLX  ( 0 . ,2.  ) +SQR*EX*CMPLX  (C,S) 

367  RAG=  (PI+ANGR)  / ( 2.0*FN) 

368  TSIN=SIN( RAG) 

369  TS=ABS ( TSI N) 

370  IF  (TS  .GT.  1 . E-5)  GO  TO  442 

371  COTA=-SQRT(  2.0  ) *F  N*S  IN(  ANGR/2 . 0 -F  N*P  I*D  N) 

372  IF(COS  (ANGR/2. 0-FN*PI*DN)  . LT  .0.0  ) COTA=-COTA 

373  GO  TO  443 

374  442  COT A=S  QRT( A) *COS ( RAG ) /T  SI N 

37  5 443  UPPI=COM*COTA*F  A 

376  DNS= (-PI+ANGR) / ( 2.0 *F N*PI) 

377  SGN=SIGN(  1.  ,DNS) 

378  N=IFIX  (ABS  (D  NS  ) +0 . 5 ) 

379  DN=SGN*N 

3 80  A=ABS  (1.0+COS(ANGR-2.0*FN*PI*DN)  ) 

381  BOTL  = 2.0*SQRT (ABS (R*A) ) 

382  EX=CEX  P(  CMP  LX  (0.0,  TP  I*  R*A  ) ) 

383  CALL  FRNELS  (C,S,  BOTL) 

384  C=SQRT (PI/2.0 )*( 0.5-C) 

385  S=  SQRT (P 1/2.0 ) * ( S- 0.5 ) 

386  F A=CMPLX  ( 0 . ,2.  ) *SQR*EX*CMPLX  (C,S) 

387  RAG=  ( P I-A  NG  R)/(  2.0*FN) 

388  TSIN=S IN (RAG) 

389  TS=ABS (TSIN) 

390  IF  (TS.  GT  .1.  E-5)  GO  TO  542 

391  COTA=  SQRT  (2.0)*FN*SIN(  ANGR/2 . 0 -F  N*P  I*D  N) 

392  IF(COS  (ANGR/2. 0 -F  N*P  I*D  N)  . LT  . 0.0  ) COTA=-COTA 

393  GO  TO  123 

394  542  COT  A=  SQRT  (A)  *COS  (RAG)  /T  SIN 

395  123  UNPI=COM*COTA*F  A 

396  DIR=UPPI+UNPI 


40 


397  RETURN 

398  END 

399  SUBROUTINE  FRNELS  (C  ,S  ,XS) 

400  C THIS  IS  THE  FRESNEL  INTEGRAL  SUBROUTINE  WHERE  THE  INTEGRAL  IS  FROM 

401  C O 0 TO  XS,  THE  INTEGRAND  IS  EXP (- J*P 1/2 . *U»U)  , AND  THE  OUTPUT  IS 

402  C C(XS)-J*S ( XS) . 

403  DIMENSION  A (1  2)  , B(  12  ) , CC  ( 12  ) , D ( L2  ) 

404  DATA  A/1. 5 95 7 69  140,  - 0.0  00  00  17  02,  -6. 8 08  56  88  54,  -0.0  00  57  63  61, 6. 9 206 91 

4 05  *90  2,-  0.0  168  98  65  7,-  3.050485660,-0.0757524  19,0 .8  5066  37  81,  - 0. 0 25  63  90  4 

406  *1,-0.150230960,0.034404779/ 

407  DATA  B/- 0.0  00  00  00  33, 4. 2 55  38  75  24, - 0.0  00  09  28  10, - 7. 7 80  02  04  00, -0.0  09  52 

408  *0895,5.075161298,-0.138341947,-1.363729124,-0.403349276,0.70222201 

409  *6,-0.216195929,0.019547031/ 

410  DATA  CC/0. ,-0.024933975,0 .000003936,0.005770956,0.000689892,-0.009 

411  *497136,0.011948809,-0.006748873,0.000246420,0.002102967,-0.0012179 

412  *30,0.000233939/ 

413  DATA  D/0. 1 99  47  11  40, 0.0  00  00  00  23, -0.0  09  35  13  41, 0.0  000  23  00  6, 0.0  04  85  14  6 

414  *6,0.001903218,-0.0  17122914,0.029064067,-  0.027928955,0.0  16497308,-0 

415  *.005598515,0.000838386/ 

416  DATA  PI/3.14159265/ 

417  IF(XS.LE.O.O)  GO  TO  414 

418  X=XS 

419  X = PI«X*X/2.0 

420  FR=0.0 

421  F 1=0 . 0 

422  K=13 

423  IF (X-4 .0 ) 10,40,40 

424  10  Y=X/4 . 0 

425  20  K=K— 1 

426  FR= ( FR+ A ( K) )»Y 

427  FI=(FI+B(K))»Y 

428  IF ( K-2)  30,30,20 
4 29  30  FR=FR+A  ( 1) 

430  F1=F I+B( 1) 

431  C-(FR*COS(X)+FI*SIN(X ))«SQRT( Y) 

4 32  S=(FR«S  IN  <X)-F  I*COS{X)  )«SQRTt  Y> 

433  RETURN 

434  40  Y=4.0/X 

435  50  K=K-1 

436  FR= (FR*CC(K) )*Y 

437  FI=(FI+D(K))*Y 

438  I F ( K-  2 ) 60,60,50 

439  60  FR=FR+CC(1) 

440  FI=FI+D(1) 

441  C=0 .5+ (FR*COS( X)+FI»S IN(X) )«SQRT( Y) 

44  2 S=0 . 5 + ( FR*S  IN(  X)  -F  I*COS(  X)  ) *S  QRT  (Y) 

443  RETURN 

444  414  C=-0.0 

445  S=- 0.0 

446  RETURN 

447  END 
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